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Diquark interaction and gaps for color superconductivity 

E. Gubankova a * 

a Institute of Theoretical and Experimental Physics, 
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Using flow equations, we derive an effective quark-quark interaction and obtain the coupled set of gap equations 
for the condensates of the CFL phase of massless Nf = 3 dense QCD. The formalism developed here enables one 
to consider more general case of nonzero s-quark mass. 



We apply flow equations to the Coulomb gauge 
QCD Hamiltonian with Nf = 3 at nonzero quark 
density. Coulomb gauge QCD Hamiltonian, V • 
A = 0, at nonzero quark density 



H 



H + H Vl 



H, 



dyn 



(i) 



includes the free Hamiltonian Ho, the instanta- 
neous interaction Hi nst describing static proper- 
ties, and the dynamical interaction Hd yn involv- 
ing gluon propagation. The free Hamiltonian is 
given by 



Ho 



dxtp(x) (—ij ■ V — fijo + m) i/j(x) 



Tr J dx (U 2 (x) + B 2 A {x)) 



(2) 



The physical degrees of freedom are the trans- 
verse gluon field A, its conjugate momentum n, 
and the quark field in the Coulomb gauge; B A is 
the abelian component of the magnetic field. The 
instantaneous interaction is given by 



Hi r} 



dxdyp a (x)V* b st (x,y)p b (y) ,(3) 



where the quark charge density is p a (x) = 
ijj(x)joT a ip(x), and in the leading order the 
kernel is diagonal ~ S ab Coulomb potential 
Vi ns t(r) = —a s /r in coordinate space or 
Vi ns t(q) = —g 2 /q 2 in momentum space. In 
high density quark matter the non-abelian con- 
tributions from three- and four-gluon interactions 
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are suppressed. Therefore we consider only the 



quark-gluon interaction, Hd yr , 



V m 



Va, 



g / dxip{x)~{ ■ A(x)ijj(x) 



(4) 



in dense QCD. 

We eliminate the quark-gluon coupling 
by the first order flow equation dV qg /dl = 
[[Ho,V qg ], H }- As a result, in the second order 
we generate the dynamical interaction between 
quarks dV dyn /dl = [[Ho,V qg ],V qg ] two— particles ; 
and the quark self-energy dT,d yn /dl = 
[[Ho,V q g],V q g] on e-particle in the channel of qq- 
pair creation/annihilation. The resulting di- 
quark interaction and self-energy include dy- 
namical terms generated by flow equations, 
Vdj/n(A — > 0) and S(A), respectively, and the 
instantaneous terms in corresponding chan- 
nels, i.e. V qq = V dyn + V inst and T, qq = 
E d?/n (A) + S mst (A) + CT(A). Here A is the 
UV cutoff connected to the flow parameter I as 
I = 1/A 2 , and CT(A) are the second order coun- 
terterms. Having the effective interaction be- 
tween quarks V qq , we allow for diquark condensa- 
tion, which can be parametrized as p| A^ 7 (p) = 

3 (± [A 8 ( P ) + |A!(p)] + iAi(p)a*a£), 

where Ai and A 8 are the eigenvalues of A^ 7 
in the CFL basis. Here i,j are flavor and a, 7 
color indices, which are substituted by p index 
in the CFL basis, 6(fc)<, = E p A£,6(fc)"/V2, A p 
are the Gell-Mann matrices for p = 1,...,8 and 
A 9 = \/2/3. Combining all the terms the high 
density effective Hamiltonian is given by H 



H 



eff 



= Ho 



V, 



qq ■ 



(5) 



2 



Using Weyl spinors, the individual terms are writ- 
ten in the CFL basis as 

#o = ^2\k-fx\ bl(k)b p (k)+J2 k a*(k)a(k) 
k, P k 

= \ J2 A p (p)e-^P\(p)b p (-p) + c.c. 



V a , 



P,P 

k.p.p.p' 



' e iHk) e - 



i<t>(P) 



x b p (k)b p (-k)b p ,(p)b p ,(-p): 



(6) 



where J^k = I dk/(2ir) 3 , and c pp ' — 1 for 
(p,p') = (1,8) or (8,1), cPP' = -2 for (p,p') = 
(8,8) and c pp ' = for (p,p') = (1,1). We 
did not include antiparticles since their contribu- 
tion is suppressed near the Fermi surface. Per- 
forming the Bogoliubov-Valatin transformation 
for quark fields from (b,w) to (y, y^), we absorb 
the self-energy into the new free Hamiltonian , 



S gg — > Hq, 



yl(k)y P (k) 

k,p 

uj{k) a\k)a(k), 



k 



(7) 



where the effective quark energy is given by 
E p (k) = ^J(k — n) 2 + Ap(fc) 2 , and the gluon en- 
ergy id includes polarization effects of a gluon 
propagating in the quark medium. The resulting 
high density effective Hamiltonian is given by 



H, 



H + V Ql 



(8) 



where V qq Eq. (g) is also subject to the BV trans- 
formation. H e ff describes dynamics of quarks 
near the Fermi surface in dense QCD. In order 
to make calculations selfconsistent we use Hq in- 
stead of H a in the flow equations. 

Further we specify the terms in the high density 
effective Hamiltonian H e ff Eq. (||) or Eq. (||). 
The kernel of the effective diquark interaction Eq. 
® is given by @ 



Vtf(k,p) = -V d p y p n (k,p) 




(9) 



where 



V d P y P Jk, P ) 



V mst (k,p) = 



2.9 2 



1 



3 6EPP'{k,p) 2 +oj M (k-p) 2 
2g 2 1 



3 u E {k-p) 2 



(10) 



and the energy difference of in- and out-going 
quarks is SEPP ' (k,p) = EP(k) - EP{p). Gluon 
energy, w(q) 2 = q 2 + M(q) 2 1 contains the De- 
bye screening mass for electric gluon, M E ~ 
g 2 p?Nf, and the Landau damping mass for mag- 
netic gluon, M M (q) 2 ~ g 2 ^ 2 N f E(q)/\q\. In the 
magnetic interaction Vd yn both factors SE 2 and 
Landau damping mass are dynamical, i.e. 

they depend on energies/momenta. Dispersion 
relation is |q| ~ E(q) when SE 2 is taken into ac- 
count, and it changes to |q| ~ E(q) 1 / 3 for Landau 
damping. Effective gluon dispersion determines 
the numerical factor in exponent for the gap solu- 
tion. Dynamical magnetic interaction generated 
by flow equations has the form — l/(q 2 + SE 2 ) as 
compared with the equal time perturbation the- 
ory result (or interaction obtained in the second 
order via Frolich transformation) — l/(q 2 — SE 2 ) 
where the energy difference has opposite sign. 
The latter interaction has an unphysical pole. To 
avoid it SE 2 is usually neglected near the Fermi 
surface. This is a valid argument when only the 
BCS singularity through the anomalous quark 
propagator is present, and the gluon propagator 
l/q 2 is replaced by a point-like interaction in the 
gap equation. However in our case of a regular in- 
teraction it is safe to keep SE 2 , which might play 
a role of the IR regulator in the collinear limit 
q~0. 

The system of gap equations for the conden- 
sates defined in Eq. (^) is given by 

Ai(p) = 8G 81 (p) 

A 8 (p) = G 18 (p)-2G 88 (p), (11) 
where the integral in the right hand side 

GPP ' M = I [J^l^M 
{P) 4 J (2^2B p (fc) 

x VPP {k,p)RPP (k,p; A) + CT(A) , 

includes the anomalous 3-d quark propagator, 
Ap/2E p , the effective gluon propagator in the 



(12) 
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form of the effective diquark interaction, V qq , 
given in Eq. (|9j), and the UV regulating func- 
tion RPp'(k,p;A) = exp(-[fc-p] 2 /A 2 ). Reg- 
ulating function generated by flow equations is 
exp (-[6EPP'(k,p) 2 +Lo M (k -p) 2 ]/A 2 ), which 

we approximate by exp (— [k — p] 2 /A 2 ) for high 
densities. It is important that the UV regu- 
lator R p P arises automatically from flow equa- 
tions that supports an underlying concept of flow 
equations being a renormalization group method. 
CT(A) is the second order local counterterm cho- 
sen from the requirement that the condesates do 
not depend on the UV cutoff, dAi/dA = for 
i = 1,8. Calculations by flow equations are done 
in the gapped theory that is reflected by the en- 
ergies E(k) containing condensates. 

System of gap equations, Eqs. ( |Tl|JT^ ), can 
be solved only numerically. However, an ap- 
proximate analytical analyses can be done by 
converting integral gap equation into differen- 
tial one [pl. Approximate solution is given by 
Aj = Aq sm(gx) with IS.^ 1 = 2&j/iexp(— Tt/2g) 
for i = 1,8, where x is a new energy variable, g 
is an effective coupling proportional to g and de- 
pending on effective gluon dispersion, and bi are 
numerical factors (&i = — 26g for color and fla- 
vor antisymmetric condensate) . Integral G in the 
right hand side of the gap equation, Eq. (p"2|), 
has the collinear (through the gluon propagator) 
and the BCS (through the quark propagator) IR 
singularities. Together these singularities give 
double logarithm structure for G and character- 
istic enhanced solution A ~ exp(— c/g) for color 
superconductivity instead of the BCS solution 
~ exp(— c/g 2 ). BSC singularity is regulated by 
the condensate A p in the energy E p . In the case of 
all three quarks are massless, the collinear singu- 
larity is cutoff by 5E(k,p) 2 for the momenta close 
to the Fermi surface |fc — n\ < A and by the Lan- 
dau damping factor Mm(&,p) 2 for the momenta 
away from the Fermi surface A < \k — fi\ < /i. 
Landau damping of a magnetic gluon is calcu- 
lated in a quark-gluon plasma, which is a nor- 
mal phase, since for large momenta gluon dis- 
tingiushes individual quarks in a Cooper pair and 
effectively 'sees' quark plasma rather than a su- 
perconducting state. For massless quarks and 



zero temperatures the characteristic double log- 
aritm in G which leads to enhanced solution is 
saturated by momenta away from the Fermi sur- 
face, and hence as first realized by Son Landau 
damping determines the gap ||. Changing pa- 
rameters may change numerical factors in the so- 
lution for the gap. In the case of massive s-quark, 
m u , md <C m s , the factor SE 2 might regulate the 
collinear singularity contributing to the conden- 
sate A. 

Generalization to nonzero temperature is done 
by inserting factor tanh(E p (k)/2T) into G, Eq. 
(p^). We assume that the effective gluon mass is 
not affected by small temperatures, since M ~ gfj, 
and T < A <C gfJ, otherwise superconducting 
condensate will melt. For the momenta away 
from the Fermi surface, A < \k — fi\ < fi, 
ta,nh(E (k) / 2T) ~ 1, hence there are no temper- 
ature effects for large energies. For the momenta 
close to the Fermi surface, \k — /i| < A, temper- 
ature distribution factor is important. Thus only 
region around the Fermi surface is affected by the 
temperature. However, this region of momenta 
does not contribute to the double logarithm of 
G and hence to the enhanced gap, but rather to 
the BCS solution (momenta close to the Fermi 
surface give single logarithm). Therefore, tem- 
perature dependence of the condensate is of the 
BCS type, i.e. A(T) ~ A(T = 0) (1 - T/T c ) 1/2 
with the connection between zero temperature 
condensate and the critical temperature given by 
T c = 0.577A(T = 0). 

In conclusion, we obtained an effective high 
density Hamiltonian for superconductivity which 
includes an effective interaction between quarks. 
Also we obtained the coupled set of gap equations 
for the condensates of the CFL phase of massless 
Nf = 3 dense QCD. The formalism used enables 
one to consider a more general case of nonzero 
s-quark mass. 
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